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Abstract. Primary and secondary masses of heavy reaction products have been deduced from kinemat- 
ics and E-ToF measurements, respectively, for the direct and reverse collisions of 93 Nb and 116 Sn at 25 
AMeV. Light charged particles have also been measured in coincidence with the heavy fragments. Direct 
experimental evidence of the correlation of energy-sharing with net mass transfer has been found using 
the information from both the heavy fragments and the light charged particles. The ratio of Hydrogen and 
Helium multiplicities points to a further correlation of angular momentum sharing with net mass transfer. 

PACS. 25.70.Lm Strongly damped collisions, 25.70.Pq Multifragment emission and correlations 



1 Introduction 

It is now experimentally established that binary dissi- 
pative processes dominate the heavy-ion reaction cross- 
section also at intermediate bombarding energies (the so- 
called "Fermi" energies) up to about 50 AMeV 
|^,[|. This evidence, based on the results of experiments 
employing very different detectors and analysis methods, 
has stirred renewed interest on the subject of dissipative 
phenomena and calls for a deeper comprehension of the 
microscopic mechanisms capable of converting such large 
amounts of kinetic energy into excitation energy. The in- 
terest is increased by the fact that at intermediate en- 
ergies dissipative processes, resembling in several aspects 
those known at low energies, coexist with new phenomena 
strictly related to the dynamics of the collision, so that 
the production of exotic and far from equilibrium pieces 
of nuclear matter may be expected. 

In this context, strongly debated topics are, among 
others, pre-equilibrium and/or mid-rapidity emission of 
light particles, the possible formation and rupture of a hot 
neck as a source of intermediate mass fragments (IMF), 
the maximum amount of thermal energy which can be de- 
posited in nuclei, the transfer of angular momentum from 
the entrance channel in the rotational modes of the frag- 
ments, the possible population of "doorway states" for 
nuclear vaporization. Another interesting and open field 
of investigation concerns the energy and angular momen- 
tum partition between the reaction partners of dissipative 
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collisions. This issue has been extensively debated in re- 
cent investigations ||[|.[7|||,|9j], but it still remains to be 
convincingly clarified. 

In previous works [||[ic],|ll| we investigated heavy- ion 
reactions at bombarding energies from 10 to 24 AMeV, 
with special effort devoted to put into evidence, in a model- 
independent way, non-equilibrium effects in dissipative col- 
lisions. In a first paper , we studied the degree of equili- 
bration between the two reaction partners at the end of the 
interaction phase employing the sequential fission chan- 
nel. We measured fission probabilities of the fragments 
produced in an asymmetric reaction ( 120 Sn + 100 Mo at 
19.1 AMeV) where a given primary mass A corresponds 
to different net mass transfers for projectile- and target- 
like fragments (PLF and TLF). The main result was that 
the curves of fission probability Pfi ss vs. fissioning mass A 
for PLF and TLF do not coincide. For a given A inter- 
mediate between target and projectile, Pfi ss for the TLF 
(which gained mass) was significantly larger than for the 
PLF, even at large TKEL (Total Kinetic Energy Loss). 
The observed effect is a clear signature of the lack of an 
overall equilibrium at the end of the interaction. 

More recently, we refined a method previously sug- 
gested by other authors fH) , thus obtaining a model in- 
dependent information on the partition of dissipated en- 
ergy J9|. The most interesting result of our work (con- 
cerning the collision 100 Mo + 120 Sn at 14.1 AMeV), was 
the observation of a sharp correlation between the num- 
ber of evaporated nucleons and the net exchanged mass 
for PLF and TLF, which we interpreted as an evidence 
for non-equilibrium excitation energy partition between 
the reaction products. This kind of correlation had been 
found and widely debated in previous works ||l^, [l3| , ^ , [l^ , 
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|l6| ] but clear conclusions on its strength and implications 
have not yet been drawn. 

Here we want to briefly summarize the terms of the 
debate on this subject. At rather low bombarding ener- 
gies (8-15 AMeV) g,0,[O|||, several experimental find- 
ings (concerning mass and charge drift, mass and charge 
variances, excitation energies of reaction products) can be 
qualitatively understood in the frame of the nucleon ex- 
change model (NEM) |2^,|2l|. In some cases a quantita- 
tive agreement can be obtained, although the experimen- 
tal variances tend to be systematically underestimated. 
The NEM assumes that energy dissipation takes place via 
a number of stochastic exchanges of single nucleons be- 
tween the two fragments when they are in contact along 
the trajectory. At the beginning of the interaction, the 
NEM predicts that the cold nuclei gain almost the same 
amount of energy due to the random nature of the ex- 
changes; as a consequence a thermal non-equilibrium de- 
velops which increases with increasing mass- asymmetry 
of the interacting nuclei. This thermal disparity tends to 
be smoothly reduced by the following nucleon transfers 
which force the dinuclear system to recover thermal equi- 
librium. As a consequence, if the duration of the contact 
is long enough (as it happens at lower impact parameters 
and at lower bombarding energies), the statistical decay 
of the two outgoing fragments should be consistent with 
the same value of the temperature (which corresponds to 
an energy partition proportional to the primary mass of 
the final fragments). 

Indeed, such a behavior has been found in many works 
|,0|l|,||,0|l|^|,^,|||4l||,|| which were aimed at 
studying the average excitation energy sharing in heavy- 
ion collisions as a function of the dissipation, commonly 
estimated by the total kinetic energy loss (TKEL); the 
word "average" means here that the data are integrated, 
for each TKEL bin, over the mass distribution of the pri- 
mary excited fragments. In these works, with increasing 
dissipation, a trend towards equilibrium partition (i.e., 
excitation energy shared in proportion to the mass of the 
fragments) is observed; however, this condition seems to 
be never reached @,||,|||6j. 



More refined experiments ||,|^,[13|,|lJ,|15j have claimed 
that the excitation energy division is correlated with the 
net mass transfer, with an excess of excitation being de- 
posited in the fragment which gains nucleons. By itself, 
the existence of such a correlation is compatible with the 
NEM framework. In fact, especially for trajectories lead- 
ing to rather peripheral collisions, it is likely that most of 
the excitation produced in each exchange comes from the 
damping of the fairly high relative velocity of the trans- 
ferred nucleon in the almost cold receptor nucleus. Thus, 
reaction products have necessarily more excitation energy 
when they have experienced mass gain than in the oppo- 
site case when they have finally lost nucleons. However, 
the observation that this correlation seems to be largely 
independent of the degree of inelasticity p|,|l5|] is another 
feature difficult to understand within the present version 
of the stochastic nucleon exchange model and deserves 
new investigation. In order to clarify this aspect of dis- 



sipative collisions we have studied the collision 93 Nb + 
116 Sn at 24.9 AMeV. 

The paper is structured as follows. In sect. || a descrip- 
tion of the detector array and of the experimental methods 
is given. Section [3] shows the experimental results which 
are divided in two parts. The first one concerns data on to- 
tal mass evaporated from the PLF as a function of its pri- 
mary mass A and TKEL (sect. 3.1); the second part deals 
with light charged particles emitted from PLF and dir ectl y 
detected with an array of plastic scintillators (sect. |3.2[ ). 
Particular care is devoted to describe the possible biases 
introduced by the measurement and analysis methods and 
the way in which they have been taken into account and 
corrected for employing Monte Carlo calculations. Details 
on this latter subject can be found in the Appendix. Fi- 
nally, sect. U presents a discussion of the results, partly 
based on a comparison with statistical model calculations 
performed with the code GEMINI §7|. 



2 Experimental set-up and methods 
2.1 The measurement 

Beams of 93 Nb and 116 Sn at 24.9 AMeV were delivered by 
the GANIL accelerator, with an excellent time resolution 
of about 550 ps (FWHM) for Nb and 350 ps (FWHM) for 
Sn. They were used to study the slightly asymmetric sys- 
tem 116 Sn + 93 Nb, both in direct and reverse kinematics, a 
method already applied in a previous experiment at lower 
bombarding energy Q. The moderate asymmetry of the 
entrance channel was necessary in order to ensure a com- 
mon range of masses for PLF and TLF even at moderate 
TKEL. 

The choice of elements Nb and Sn was motivated by 
beam and target feasibility reasons, with the additional 
constraint of studying a system not too far from 100 Mo 
+ 120 Sn, previously studied at 14 AMeV at GSI §. The 
particular isotope of Sn was thus a compromise between 
the technical need for a sufficiently large natural abun- 
dance (to make it usable as a beam source) and the re- 
quirement of having an N/Z ratio (1-32) similar to that of 
93 Nb (1-27), in order to reduce the possible role of isospin 
equilibration. Both targets consisted of foils of 93 Nb and 
isotopically enriched 116 Sn, with a thickness of about 200 
/xg/cm 2 . 

During most runs, the target was in a tilted position, 
with the beam impinging on it at an angle of about 45°. 
This choice, while not appreciably affecting the measure- 
ment of the fast forward-going PLF, was of utmost im- 
portance for a good measurement of the coincident TLF, 
especially for the slow ones emitted close to 90° (as in case 
of events with moderate energy dissipation) . In this way 
one strongly reduces the average thickness of target ma- 
terial passed through by the TLF, drastically decreasing 
the perturbation of its velocity vector due to energy- and 
angle-straggling. 

Shorter runs with other targets were also performed 
at the same bombarding energy. In particular, as it will 
be explained in the following, the data collected for the 
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two symmetric systems 116 Sn+ 116 Sn and 93 Nb+ 93 Nb were 
used for checks and corrections. Further information was 
obtained from short runs with the strongly asymmetric 
systems 116 Sn+ 58 Ni and 116 Sn+ 197 Au at the same energy. 

Finally, for calibration purposes, data were recorded 
also for the system 93 Nb+ 93 Nb at a bombarding energy as 
low as 3.86 AMeV, obtained by switching off the second 
cyclotron of GANIL. In this way the whole set-up was 
illuminated with elastic scattering events. 



2.2 The set-up 

The experiment is based on the determination, as a func- 
tion of TKEL, of the total mass evaporated from excited 
PLF, as given by the difference between their primary 
mass (obtained, via the kinematic coincidence method, 
from the velocity vectors measured by gas detectors) and 
secondary mass (via additional measurement of the kinetic 
energy by means of Silicon detectors) . Additional informa- 
tion on the dissipative collisions was also obtained from 
the direct measurement of the light charged particles by 
means of scintillation detectors. 



2.2.1 The gas detectors 

An array of 12 large-area position-sensitive parallel-plate 
avalanche detectors (PPAD) was mounted inside the Nau- 
tilus scattering chamber. The detectors were arranged on 
three planes of four detectors each, in an almost axially 
symmetric geometry around the beam direction (see fig. 1 
of ref . Q ) . With respect to previous experiments at lower 
bombarding energies JD>[l9], the most forward plane was 
moved at a larger distancefrom the target (about 250 cm), 
in order to increase the flight-path, and the resulting larger 
dead-region between the forward and the middle plane was 
partially reduced by shifting two of the middle plane de- 
tectors toward the beam axis. Thus, in this experiment, 
the 12 gas detectors (each with an active area of 300x300 
mm 2 ) covered about 65% of the forward hemisphere. 

The PPAD detected heavy (Z^IO) reaction products, 
with about 100% intrinsic efficiency. They measured both 
the position of impact and the time-of-ffight (with re- 
spect to the bunched beam) of the reaction products, thus 
yielding their velocity vectors. The position resolution of 
the PPAD was about 3.5 mm (FWHM) and the overall 
timc-of-ffight resolution (including the contribution of the 
beam) about 750 and 600 ps (FWHM) for the runs with 
Nb and Sn beam, respectively. 

In the experiment, the most critical operating condi- 
tions were those of the most forward PPAD, due to the 
high counting rate of elastic and quasi-elastic products. 
Thus, typical beam currents of 0.1-0.5 nA (with charge 
states 31 + and 37 + for the two beams) were used to limit 
the maximum counting rate of the forward detectors to 
about 15-20 • 10 3 counts/s. Moreover, in order to reduce 
the load on all detectors due to the flux of electrons ex- 
tracted by the beam while passing through the target, a 



positive high voltage of 40-42 kV was applied to the elec- 
trically insulated target holder. More details on the gas 
detectors can be found in [Q. 

2.2.2 The silicon detectors 

Two identical arrays of 23 ion-implanted Silicon detectors 
each were mounted behind two of the four most forward 
PPAD. The Silicon detectors covered laboratory polar an- 
gles ranging from about 2° to about 7°, so that most of 
them were located below or around the grazing angles for 
the studied reactions, which vary from about 3.5° for the 
116 Sn+ 58 Ni system to about 8.6° for the 116 Sn+ 197 Au one. 
The first 18 Silicon detectors of each array (manufactured 
by Eurisys Mesures) had an active area of about 30x30 
mm 2 and a thickness of about 500/im, sufficient to fully 
stop the quasi-elastic PLF at these bombarding energies. 
They covered the region at smaller angles, approximately 
between 2° and 5.5°. In each array, larger angles up to 
about 7° were covered by 5 Silicon detectors (purchased 
from Micron Semiconductors for previous experiments) 
with an active area of 50x50 mm 2 and a thickness of 
300/im. In order to fully stop quasi-elastic PLF, they were 
mounted in a tilted position with respect to the direction 
of the incoming particles, so that their effective thickness 
was increased to about 420/im and their effective area cor- 
respondingly reduced to about 35 x 50 mm 2 . 

2.2.3 The scintillation detectors 

For the measurement of the light charged particles (LCP) 
we took advantage of the scintillator array "Le Mur", 
mounted at small angles on the closing cup of the NAU- 
TILUS scattering chamber, behind the gas and Silicon de- 
tectors. This device consists of 96 pads of fast plastic scin- 
tillator NE102, 2 mm thick, mounted in 7 circular rings 
centered on the beam axis, with a threshold of about 3.2 
AMeV for protons and a-particles. A detailed description 
of the geometry and performance of this device can be 
found in pq ], "Le Mur" allows a clean Z-identification 
of fast light reaction products punching through the thin 
scintillator material. For protons and a-particles this hap- 
pens for energies greater than sal3.5 AMeV. 

Due to the primary need to optimize the time of flight 
resolution, the target holder had to be installed at one 
extreme of the scattering chamber, thus increasing the 
flight-path but significantly reducing the angular accep- 
tance of the wall (polar angles from about 2° to about 
18.5 °). Moreover, the presence of our apparatus between 
target and "Le Mur" (with absorbing materials like de- 
tector frames, cables and supporting structure) produced 
large and complicated shadows on several plastic pads. In 
the analysis of the data it was then necessary to select a 
subset of pads which were reasonably clean and free from 
shadows. A first selection then was operated by cutting 
away all the pads in which distortions or anomalous lack 
of yield in the ridges for Z=l and Z=2 particles were a 
clear evidence of shadowing. Some detectors of the inner 
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rings which were contaminated by the scattering of beam 
particles on the last collimator, as well as detectors in 
which poor resolution did not allow a clean separation of 
Z=l and Z=2 particles, were discarded too. Finally, only 
33 out of the original 96 pads were retained in the analy- 
sis. In the identification matrices also fast light fragments, 
such as Lithium or Beryllium, were visible above the ridges 
of 7j—1 and Z=2 particles, however the cumulated statis- 
tics was low and they will not be considered in the present 
work. 

The geometric acceptance of "Le Mur" was such that, 
for our slightly asymmetric systems, most of the detected 
LCP originate from PLF decay, except for very central 
collisions. We even strengthened this geometric selection 
by cutting away in the analysis all slow particles stopped 
in the pads, the amount of rejected particles being about 
15-20% of the total. 



2.3 Calibration and correction procedures 

2.3.1 Kinematic coincidence method 

A refined version of the kinematic coincidence method 
(KCM) (29) was applied and from the velocity vectors 
measured by the PPAD primary (pre-evaporative) quan- 
tities -in particular the masses of the heavy fragments 
and TKEL- were deduced event-by-event. Although this 
analysis can be applied to events with 2, 3 or 4 heavy frag- 
ments in the exit channel, the main interest of this work 
is concentrated on the binary channel. Therefore, unless 
otherwise explicitly stated, in the following we will refer 
only to measured binary events in which two heavy frag- 
ments (with charge Z>10) were detected by two PPAD 
and one of them was also stopped in one Silicon detector. 

This version of the kinematic coincidence method al- 
lows to exploit the two-fold redundancy of the available 
experimental information for 2-body events. Indeed, it is 
based on the minimization of J2i=i 2 m i l( v i XP — v «)| 2 
under constraint of mass and momentum conservation. It 
gives not only the best estimates of the unknown primary 
masses m, , but also optimized "improved values" v< of the 
measured velocity vectors v| xp , such that the conservation 
laws are exactly satisfied event-by-event. On the basis of 
statistical arguments |^9| , for binary events the distribu- 
tion of the minimized quantity Zikcm is expected to be ap- 
proximately shaped like a x-distribution with 2 degrees of 
freedom. The width of the distribution increases with in- 
creasing perturbation of the 2-body kinematics (i.e., with 
increasing dissipation) . 

2.3.2 Correction for misalignment of the set-up 

A precise knowledge of the position of the detectors with 
respect to beam direction is mandatory for a kinematic re- 
construction, especially in case of strongly forward-peaked 
kinematics. Therefore special care was devoted to the mea- 
surement of the geometry of the set-up and to its align- 
ment along the nominal beam direction. The geometry of 



the detectors was measured optically by means of a sur- 
veying instrument (Total Station) mounted in the target 
position. The nominal accuracy of the instrument for an- 
gle and distance measurements (the last one obtained via 
an infrared distantiometer and corner cubes on the detec- 
tors) is about 5" and ±2 mm, respectively. Due to the 
uncertainty (±0.5 mm) in the positioning of the instru- 
ment with respect to the target center and of the corner 
cubes with respect to the detectors, angular accuracies of 
the order of 6' for the detectors closer to the target and of 
40" for the more distant ones are estimated. The absolute 
distance accuracy is the nominal one of ±2 mm. 

The beam was neatly focused onto the target center 
with the help of an alumina (AI2O3) plate mounted on 
the target ladder: the beam was steered until its lumines- 
cent spot almost disappeared into a hole (of 2 mm radius) 
drilled in the center of the alumina. The coincidence of 
actual direction of the beam with the one optically deter- 
mined during the assembly of the setup was checked at 
run-time by a beam profiler (consisting of a small remov- 
able grid-detector) located on the beam axis about 3 m 
downstream of the target. 

In the off-line analysis, the position of the detectors 
was checked with the data of elastic scattering of 93 Nb 
projectiles on light targets ( 12 C, 27 Al and 55 Mn), taken 
at the lower bombarding energy of 3.86 AMeV. In fact, 
in asymmetric systems and reverse kinematics, the elas- 
tic scattering of heavy projectiles presents a limiting an- 
gle, which in our case was 6*}™" w 7.4°, 16.9° and 36.3° 
for the three light targets, respectively. The presence of a 
limiting angle produces a circular ridge in the cross sec- 
tion d 2 er/d#i a b d<f>, characterized by a sharp drop of the 
intensity towards larger laboratory angles, thus allowing 
an easy consistency check of the geometry of the detectors. 

Finally, the actual alignment of the beam axis with re- 
spect to the optical axis was checked in the off-line analysis 
by means of the elastically scattered projectiles which par- 
tially irradiated the four most forward PPAD in all studied 
reactions. To do so, it was assumed that the beam hits the 
target center, but its direction may present a slight devi- 
ation from the optical axis. The proper correction would 
then consist in a small rotation of the reference frame 
around the target center. 

For each run, the angular distributions of elastically 
scattered projectiles hitting the four most forward PPAD 
(in clean equal windows of the azimuthal angle <fi) were 
simultaneously fitted with the Rutherford cross section, 
0Ruth(#iab)- The "best" correction was determined by \ 2 
minimization with respect to the applied reference frame 
rotation. For all runs, such a correction amounted to 0.1° 
at most (less than 4 mm in the plane of the most forward 
detectors) . 



2.3.3 Correction for pulse height defect 

The Silicon detectors measured the kinetic energy of heavy 
fragments passing through the PPAD. From this kinetic 
energy, and from the time-of-flight (ToF) of the PPAD, 
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secondary masses for PLF were deduced employing an it- 
erative procedure which includes corrections both for the 
"pulse height defect" (PHD) in Silicon and for the energy 
lost in the PPAD and in the dead layers of the Silicon chip. 
Quantitatively, both corrections represent a small fraction 
of the total kinetic energy for most ions. As an example, 
for elastically scattered 116 Sn or 93 Nb ions at 25 AMeV, 
the total energy-loss before reaching the active region of 
the Silicon chip amounts to about 1% of their energy, and 
the correction for the PHD in the bulk of the chip is at 
most comparable. However, since good accuracy in the 
energy measurement was a strong need for our analysis 
purposes, special effort has been devoted, both during the 
measurement and in the off-line analysis, to obtain a pre- 
cise evaluation of the energy released in the Silicon detec- 
tors. Experimental values of PHD (extracted from data) 
were taken into account in the calibration and analysis. 
For a careful discussion of these aspects we refer to the 
technical paper j$o|. Here, we only recall that the over- 
all accuracy of the energy calibration is estimated to be 
around 0.5%. 

In order to judge about the quality of the measure- 
ment of both primary and secondary mass of PLF, some 
experimental results are shown in fig. [l]. Part a) presents 
the kinetic energy spectrum of PLF, measured with one of 
the Silicon detectors (located at « 2.7°) in binary events of 
the 116 Sn+ 93 Nb reaction at 25 AMeV. One can clearly see 
the elastic peak, with a resolution of 7-8 MeV (FWHM), 
which represents a typical value for all the Silicon detec- 
tors. Gating on this peak (for each Silicon detector) gives 
an easy way to select (or exclude) elastic events in the 
analysis. 



2.3.4 Corrections for position-dependent time-of-flight 

Figures |lb,c,d show the resolutions in time-of-flight, pri- 
mary and secondary mass, respectively, for clastic events 
selected by gating on the elastic peaks in the Silicon en- 
ergy spectra. Here and in the following, the time-of-flights 
measured by the PPAD have been corrected for the depen- 
dence of signal propagation time on the position of impact 
on the PPAD. This correction was indeed crucial for ob- 
taining the displayed good resolutions and was deduced, 
for all PPAD, from the analysis of pure elastic events pro- 
duced in collisions at low bombarding energy ( 93 Nb+ 93 Nb 
at 3.86 AMeV). 

Figure |l|b presents the distribution of the differences 
between the experimental time-of-flight of the PPAD and 
the one calculated via elastic scattering kinematics from 
the deflection angle measured by the PPAD. All events 
with projectiles hitting any one of the Silicon detectors 
behind one of the forward PPAD have been added. As can 
be seen, the overall time-of-flight resolution was indeed 
very good (At < 600 ps FWHM), thanks also to the high 
quality of the beam buckets delivered by the GANIL staff. 

The primary mass distribution for these elastic events, 
as obtained from the kinematic coincidence method, is 
shown in fig. ^jc. It is worth noting that the good mass res- 
olution of about 2 amu (FWHM) depends also on the good 
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Fig. 1. a): Experimental kinetic energy spectrum measured 
with one of the Silicon detectors for binary events in the 
116 Sn+ 93 Nb system at 25 AMeV. In this figure no subtrac- 
tion of the ions scattered on the grid wires of the PPAD has 
been performed. In the portion of spectrum displayed in the in- 
set, the continuous line represents all the 2-body events, while 
the dashed line is the estimated background due to the ions 
which pass through the grid wire s of th e PPAD before reach- 
ing the Silicon detector (see sect. 2.3.5| in the text). The data 
presented in the remaining panels of the figure have been gated 
by windows on the elastic peaks of the detectors, b) : Distribu- 
tion of the differences between measured and calculated time- 
of-flight for elastically scattered projectiles; all the Silicon de- 
tectors behind one of the PPAD were taken into account, c): 
Primary mass distribution of elastically scattered projectiles, 
obtained from binary events by means of the kinematic coinci- 
dence method, d): Secondary mass distribution obtained from 
the energy measured in the Silicon detectors and time-of-flight 
measured with the PPAD. 



angular resolution for the very slow target recoil accom- 
panying the fast projectile in these binary events. Finally, 
for the same sample of elastic events as in parts b) and c), 
fig. shows the secondary mass distribution, which is ob- 
tained from quantities related solely to the PLF detection. 
As expected, in this case the mass resolution is dominated 
by the time-of-flight resolution which accounts for nearly 
90% of the measured value of 4.2 amu (FWHM). 



2.3.5 Effects of the wires of the gas detectors on PLF 
energy measurement 

A tricky effect is connected with the measurement of the 
secondary mass of ions passing through the position sensi- 
tive PPAD. The grids of wires of the PPAD anodes, from 
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which the impact position is deduced, represent an in- 
homogeneous dead-layer for the impinging ions. At these 
high beam energies, energetic reaction products are not 
stopped in the wires, even not when they cross the whole 
diameter. In this case, one obtains the correct information 
on ToF and position from the PPAD, but the measured ki- 
netic energy in the Silicon detectors is strongly degraded. 
This effect concerns a small fraction of events, about 4%, 
determined by the geometric "cross section" of the two or- 
thogonal grids (each made of 20 fim diameter wires with a 
pitch of 1 mm). Due to the large intensity of the Ruther- 
ford cross section, a large amount of these kind of events 
corresponds to degraded elastically scattered projectiles, 
simulating inelastic events. They are responsible for the 
peaked distribution slightly above 2 GeV which can be 
clearly seen in fig. |l|a. 

However, not only the energy of the elastic projectiles, 
but also that of other inelastic products gets similarly de- 
graded when passing through the PPAD, thus producing 
an unknown background of bad events, which one would 
like to single out and reject. Every attempt to eliminate 
them on the basis of correlations between measured quan- 
tities (like e.g. ToF and AE from PPAD and E from the 
Silicon detectors) introduced spurious cuts as a side-effect. 
Thus we devised an analysis procedure to get an average 
estimate of such a background. All reaction products de- 
tected in the Silicon detectors (including the elastic ones) , 
were "passed" through a "simulated" grid, namely the 
point of impact on the wire was randomly chosen, the 
thickness of material calculated and the energy accord- 
ingly degraded (also the occurrence of double impacts on 
both grids was taken into account with the proper weight). 
The obtained "degraded" data were re-analyzed and the 
so estimated background was subtracted from all results 
based on the information of the Silicon detectors, with a 
suitable normalization chosen so as to remove the spurious 
peak in the energy spectra of fig. |l|a. 

The quality of the correction can be judged from the 
inset in fig. [l]a, which shows a part of the kinetic energy 
spectrum. The continuous line represents all the 2-body 
events, while the dashed line represents the background of 
nuclei slowed down by the wire grids of the PPAD, as esti- 
mated with the above described procedure. Similar results 
are obtained for other detectors and for the 93 Nb+ 116 Sn 
reaction. 



2.3.6 Calibration of the scintillators 

An example of the correlation light-output vs. time-of- 
flight is presented in fig. ||for the reaction 116 Sn + 93 Nb at 
25 AMeV. In this correlation each species of light charged 
particles gives origin to a ridge divided into two branches 
by a cusp. To the left of the cusp there are fast parti- 
cles punching through the thin scintillator, while to the 
right there are slower particles which are stopped in the 
material. The position of the cusp was used as a calibra- 
tion point for the time-of-flight scale. In case of punch- 
ing through particles, the light output is related to the 
energy lost in the material and a clean Z-identification 
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Fig. 2. Experimental scatter-plot of light-output vs. ToF 
for light charged particles measured with one of the pads of 
the scintillator array "Le Mur" for the system 93 Nb+ 116 Sn at 
25 AMeV. The particles are detected in coincidence with two 
heavy reaction products in the PPAD. The ridges of energetic 
punching-through Z=l and Z=2 particles are clearly visible. 



can be achieved, in spite of quenching effects. The two in- 
tense branches on the left of fig. || correspond to energetic 
Z=l and Z=2 particles. In case of stopped particles, the 
light output is related to their total energy, but strong 
quenching effects prevent a clean mass-identification. For 
example, it is well known that stopped a-particles become 
mixed up with stopped deuterons and tritons. 



2.3.7 Background of incompletely detected events 

As stated before, in this paper we focus on dissipative 
binary events, in which two heavy fragments were de- 
tected by the PPAD and one of them (usually the PLF) 
was also stopped in one of the Silicon chips. Similarly to 
previous works Q, also in the present case true binary 
events still represent a maj or part of the total reaction 
cross section (see sect. 3.1.1). However, a sizeable amount 
of higher multiplicity events is also produced, especially 
at high TKEL values. Therefore, because of the incom- 



plete geometric coverage of the set-up, a certain fraction 
of the detected 2-body events are not true binary events, 
but rather partially detected events of multiplicity greater 
than two (mainly 3-body events). Following a procedure 
described in detail elsewhere j2^, we estimated and sub- 
tracted this background from the data. 

The continuous histograms in fig. || show, for some 
bins of Total Kinetic Energy (TKE) in the reaction 93 Nb 
+ 116 Sn at 25 AMeV, the experimental distributions of 
^kem, the quantity minimized in the kinematic coinci- 
dence analysis. It has been verified that the histograms are 
approximately shaped like a ^-distribution with 2 degrees 
of freedom (29|, except for the long flat tails. The dot- 
ted histograms give the (estimated) 3-body background 
obtained by analyzing the measured 3-body events (with 
one fragment in one Silicon detector) as if they were binary 
events, after having randomly removed one of the two frag- 
ments which did not hit the Silicon detector. They have 
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Fig. 3. Spectra of the Zikcm for the quoted bins of TKE for 
the reaction 93 Nb+ 116 Sn at 25 AMeV; the full curve refers to 
the measured 2-body events while the dashed one represents 
the estimated background due to incompletely detected 3-body 
events. The quoted percentage of the background contribution 
refers to the region below Zikcm ~ 8 mm/ns. 



2.3.8 Correction for intrinsic and geometric efficiency 

In order to obtain meaningful results, all experimental dis- 
tributions need to be corrected for "instrumental" and 
physical effects, due not only to the finite resolution of 
the measurement and possible biases of the analysis, but 
also to the smearing of the particle evaporation process. 
As an example, the average values deduced from nonuni- 
formly distributed variables may be severely biased, due to 
finite resolution effects, unless proper corrections are ap- 
plied (see, e.g., the comment about angular distributions 
in |2j| and the correction of mass distributions in pl[). 

Moreover, it is worth noting that the quantities of 
physical interest (like, e.g., masses, angles or dissipated 
energy), although derived from truly uncorrelated param- 
eters measured by the detectors (e.g., timc-of- flight, x-y 
position, deposited energy and so on), acquire a certain 
degree of correlation pll, which must be corrected for. 

In general the corrections are rather involved and may 
be worked out analytically only in very simple cases. There- 
fore in the present work the experimental results were cor- 
rected via extensive Monte Carlo simulations, modeling 
the dissipative collision followed by an evaporative emis- 
sion and incorporating as realistically as possible the re- 
sponse of the setup, finite resolution effects and distortions 
of the analysis method. 

The parametrization of the dissipative collision was 
tuned in such a way that the simulated distributions, af- 
ter passing the experimental filter, reproduced the exper- 
imental ones. For the evaporative step, thermal emission 
(leading to Maxwellian energy distributions) was assumed, 
until the excitation energy of the emitter was exhausted. 
The multiplicity of light particles was tuned on the re- 
sults of statistical model calculations with the code GEM- 
INI [^7| , while the multiplicity of intermediate mass frag- 
ments (IMF) had to be so mewha t enhanced to reproduce 
the data (see later in sect. 3.1.3). More details about the 



simulation can be found in the Appendix and in ref. 



3 Data analysis and experimental results 



been normalized (by a single normalization factor for all 
TKEL bins) to the experimental distributions in the tail 
region. As can be seen, the shape of the long tails at high 
^kem values is well reproduced, thus lending further sup- 
port to its interpretation in terms of 3- (or more-)body 
background. 

Only events in the window A^ cm < 8 mm/ns, corre- 
sponding to the peak of the distribution, were consid- 
ered in the analysis. Moreover, the contribution of the 
estimated background was subtracted from the final spec- 
tra pjJ. This contribution grows with increasing dissipa- 
tion until, at TKEL of 600-700 MeV, it accounts for almost 
30% of the measured events. This is the reason why, in the 
following, we shall limit our analysis at 2-body events with 
TKEL?C700 MeV. 



An estimate of the sharing of excitation energy between 
the two reaction products of a binary dissipative collision 
can be performed by detecting also the emitted light par- 
ticles. An alternative, and in many respects complemen- 
tary, method consists in the simultaneous estimate of the 
primary (or pre-evaporative) mass A and secondary (or 
post-evaporative) mass A scc of the products of a binary 
collision, so that the total number of emitted nucleons 
can be obtained from the difference between these two 
values, A A = A — A scc . At the cost of no information 
whatsoever on the various particle species, this method 
yields the global number of nucleons (irrespective of en- 
ergy, emission angle and species of particle) which were 
lost by the investigated fragments. 

This method was applied in the past to measurement 
of PLF from rather asy mmetric systems studied in direct 
kinematics only |l^,|2j||l5| and required a detailed and not 
trivial comparison between the experimental results and 
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evaporation calculations. To avoid relying on model cal- 
culations (which become increasingly uncertain with in- 
creasing excitation energy), we aimed at comparing not 
the data with a model, but directly two sets of experi- 
mental data. 

With an asymmetric colliding system, one might com- 
pare the two event samples in which reaction products of a 
given mass A are PLF or TLF, this fact implying different 
"histories" (gained or lost nucleons). To overcome the se- 
vere experimental difficulties (like threshold effects, poor 
resolution, and critical dead layer corrections) which im- 
pede the measurement of the secondary mass of the TLF 
with sufficient accuracy, we devised the alternative ap- 
proach of measuring the secondary mass of the PLF only, 
however studying the same asymmetric collision both in 
direct and reverse kinematics. This approach gives also 
the additional bonus that the efficiencies for the detection 
of the PLF, being quite similar for the two kinematics, 
practically do not affect the result of the comparison. 



3.1 Experimental results on total evaporated mass 

3.1.1 General features 

The two most abundant event types are those with two or 
three heavy fragments in the exit channel. Their relative 
importance is shown in fig. |], where the two yields are 
plotted as a function of TKE for the reaction 93 Nb+ 116 Sn 
at 25 AMeV, after correction for the efficiency of the set- 
up by means of Monte Carlo (MC) simulations. As can 
be seen, the binary exit channel accounts for most part of 
the total reaction cross section, but the ternary channel is 
dominant at low TKE values. 

The gross features of the binary events resemble those 
of the well known decp-inclastic reactions at lower bom- 
barding energies (fig. |). The distribution of TKE extends 
continuously from quasi-elastic energies down to very low 
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Fig. 4. Distributions of 2- (solid line) and 3-body (dashed line) 
events as a function of TKE for the reaction 93 Nb+ 116 Sn at 
25 AMeV. The data are corrected for the different detection 
efficiencies of the two exit channels via Monte Carlo simula- 
tions. 



values, roughly corresponding to the Coulomb repulsion 
of two nuclei at contact. 

The angular distribution of PLF is strongly focused 
near or below the grazing angle for events in a wide range 
of TKE values (down to about 500 MeV) , corresponding 
to partly relaxed events. Only fully relaxed events present 
a broad and nearly flat angular distribution, which is sug- 
gestive of a possible orbiting behavior. 

Concerning the primary mass distributions, at the end 
of the interaction the two reaction products maintain, on 
the average, their original mass value, but the variance 
of their mass distribution undergoes a rapid growth when 
passing from peripheral collisions to more central ones. 
This latter fact can be seen in fig. || which shows, as a 
function of TKEL, the variance a\ (open symbols) of the 
primary (MC corrected) mass distribution of PLF from 
the reaction 93 Nb+ 116 Sn at 25 AMeV, measured in direct 
and reverse kinematics. The variances are obtained from 
Gaussian fits to the primary mass distributions obtained 
with the kinematic coincidence method. 
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Fig. 5. Experimental correlations d 2 a/d TKE d 0p™ F 
(Wilczynski-Plot, upper part) and d 2 a / AT K E AAp^ F 
(Diffusion-Plot, lower part) for the system 93 Nb + 116 Sn at 
25 AMeV, not efficiency corrected. 
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3.1.2 Mass variances, AA, A 



ken 



In particular, as a good estimate of primary masses was 
of paramount relevance in this work, a great effort was 
devoted to a better understanding of possible distortions 
which may affect the measurement of primary masses. 
Although we correct for these effects with Monte Carlo 
methods, it can be shown that the distributions of re- 
constructed primary mass are dominated by the physical 
width and not by the "instrumental" resolution, so that 
the corrections are small. For comparison, fig ^ shows the 
contribution to the variances due only to the evaporation, 
detection and reconstruction effects (stars in the figure). 
It was estimated by assuming, in the Monte Carlo simu- 
lation, that the primary masses in the exit channel were 
those of projectile and target. As can be seen, this con- 
tribution to the observed mass variance is always of the 
order of 25% or smaller. Thus, from the simple formula 
proposed in J3lj , one can get a quick idea about the errors 
affecting the uncorrected primary-mass values. Assuming 
A=93 for the most probable exit mass of the PLF, an 
uncorrected value of A=103 would be overestimated by 2 
amu, at most. 

In fig. [3] we also show (full line) the mass variances 
obtained from Monte Carlo simulations, when using the 
parameterization proposed by Gralla et al. fl32| for simi- 
lar mass systems at lower bombarding energies, with the 
values of the parameters tuned to better reproduce the ex- 
perimental widths at 25 AMeV (for the exponential slope 
parameter in eq. 1 of [32] j, a value of -4.3 fi./MeV was used 
instead of -5.27). 

In dissipative collisions, large amounts of excitation 
energy are deposited in the outgoing primary fragments, 
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Fig. 6. Experimental variances a\ of the primary mass dis- 
tributions of PLF, as a function of TKEL, for the collision 
93 Nb+ 116 Sn at 25 AMeV, both in direct and reverse kinemat- 
ics. PLF masses and TKEL were estimated with the kine- 
matic coincidence method and the variances were obtained 
with Gaussian fits after efficiency correction. The full line 
shows the parametrization of the variances used in the Monte 
Carlo simulations for the efficiency corrections. The stars rep- 
resent the contribution to the variances due to evaporation, 
detection and reconstruction, as deduced from Monte Carlo 
simulations. 
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Fig. 7. a): Distribution of the difference A A = A — A sac 
between primary and secondary mass of the PLF as a function 
of TKEL for the reaction 116 Sn + 93 Nb at 25 AMeV. The data 
are not efficiency corrected. The shape of the distribution in 
the (quasi-)elastic region is due to the finite resolution and 
to the correlations arising from the use of the measured time- 
of-flight for the estimation of masses and TKEL. b): Average 
value of AA for different windows of TKEL in the collisions 
116 Sn+ 93 Nb at 25 AMeV (full circles) and 120 Sn+ 100 Mo at 
14 AMeV (open circles). The vertical bars represent the second 
moments of the distributions. 



which de-excite mainly by evaporation of light particles. 
This feature is clearly apparent in fig. ^a which shows, 
for the reaction 116 Sn+ 93 Nb at 25 AMeV, the measured 
correlation (not Monte Carlo corrected) between the to- 
tal number of nucleons emitted by the PLF, A A = A — 
A sec , and TKEL. These results are integrated over all PLF 
masses. The weak tail extending towards large values of 
A A at small TKEL, is what remains of the elastic pro- 
jectiles p unchin g through the wires of the gas counters 
(see sect. 2.3.5) after the subtraction of this background. 
This plot looks quite similar to that reported in fig. 1 of 
the paper by Kwiatkowski et al. p5| . In the present case, 
however, the de-excitation of nuclei with much larger exci- 
tation energies can be studied, at TKEL values as large as 
«700 MeV, where up to about 50 nucleons are lost by the 
PLF. By cutting the bidimensional distribution of fig. [?]a 
into vertical slices, one obtains the average values of AA as 
a function of TKEL presented by the full circles in part b) 
of the same figure, with the vertical bars representing the 
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second moments. For comparison, the same figure shows 
as open circles also the AA of PLF produced in the reac- 
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tion 12U Sn- 



J Mo at 14 AMeV 



The data at the two 



bombarding energies practically coincide, as it has to be 
for an evaporative process if TKEL is a good estimate of 
the total excitation energy of the system and AA a good 
estimate of the excitation energy of the PLF. This fact can 
be taken as an indirect evidence that, at 25 AMeV and for 
relatively large impact parameters, pre-cquilibrium effects 
still do not play a major role. In fact, in nearly symmetric 
collisions of heavy nuclei, on average, the center-of-mass 
of the interacting system after pre-equilibrium emission is 
not expected to be appreciably different from the center- 
of-mass of the entrance channel. The main effect of pre- 
equilibrium on the kinematic reconstruction should then 
be an overestimation, by a common scaling factor, of the 
true values of E cm , TKE and A. As a consequence the 
difference AA = A — A scc should be much more strongly 
affected than TKEL= E cm - TKE and this seems not to 
be the case in fig. [?]b. 



3.1.3 Intermediate mass fragments 

Figure H shows the shape of the experimental distributions 
of A A (full histograms) for several TKEL windows. The 
spectra were obtained by projecting TKEL-slices of fig. 
onto the vertical axis, but with the additional condition 
that the primary mass A be in a narrow window around 
the mass of the projectile, namely at 116±2. With in- 
creasing TKEL, the maximum of the distributions moves 
to larger values and the width increases, while an evident 
tail towards large values of A A is present at all TKEL. 
This can be considered as an indirect evidence for the 
emission of Intermediate Mass Fragments (IMF) which 
may be important especially in the tail of the distribu- 
tions towards large AA values. With respect to light par- 
ticles, IMF have a larger binding energy per nucleon and 
in a thermal process they are emitted with lower kinetic 
energy per nucleon. Therefore, for a given total emitted 
mass, they should be less effective in removing excitation 
energy. Thus, a nucleus of a given mass and excitation en- 
ergy which evaporates one or more IMF is likely to end up 
with a smaller secondary mass (and hence a larger AA). 

Indeed, a good reproduction of the experimental dis- 
tributions was obtained when including in the evaporation 
step of the Monte Carlo simulation not only the emission 
of light particles, but also of IMF (actually we took just 
Li and C nuclei as representative of all IMF). Their mul- 
tiplicity, linearly rising with excitation energy, was tuned 
to reproduce the data and the final results are shown by 
the dashed histograms in fig. ||. The adopted values corre- 
spond to total multiplicities of IMF about 3 times larger 
than predicted by the statistical code GEMINI. In order 
to put into evidence the importance of IMF for reproduc- 
ing the tail on the right side of the distributions, we run 
also a simulation featuring the statistical emission of light 
particles only (dotted histograms). While the position of 
the most probable value and the shape of the left side of 
the experimental distributions are rather well reproduced, 
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Fig. 8. Shape of the distributions of AA for different bins of 
TKEL from the reaction 116 Sn + 93 Nb at 25 AMeV. Shown 
are the spectra for the experimental data (full histograms), the 
results of a Monte Carlo simulation featuring only evaporation 
of light particles (dotted histograms) and the results of a Monte 
Carlo simulation including also the emission of IMF, tuned to 
reproduce the experimental data (dashed histograms). Note 
the changes of abscissa. 



the tail on the right side is not reproduced at all. At low 
TKEL it underestimates the experimental data and with 
increasing TKEL the obtained distributions tend to be- 
come even more symmetric. 

Of course mechanisms other than statistical emission 
may contribute and indeed, at somewhat larger bombard- 
ing energies, IMF emission from the "neck" region has 
been observed 



34 1). However, the effect 
le AA distribution is 



(see, e.g., refs. [|33 
of IMF emission on the shape of t 
indirect and not very sensitive to the assumed mechanism, 
so no attempt to refine the simulations was done. 

Indirect evidence for an increased emission of IMF can 
be obtained not only from the shape of the total emitted 
mass AA, but also from the perturbation of the 2-body 
kinematics. In fact, the recoil effects due to the IMF emis- 
sion represent a sizeable perturbation and are visible also 
in t he inc reased width of the distribution of A^ cm (see 
sect. 2.3.1). Figure H shows the shape of the experimental 
distribution of Z\kcm (full histogram) in one bin of TKEL, 
compared with the result of the Monte Carlo simulation 
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A kcm (mm/ns) 

Fig. 9. Shape of the distribution of zikcm in one bin of TKEL 
from the reaction 116 Sn + 93 Nb at 25 AMeV. As in the previ- 
ous figure, full, dotted and dashed histograms represent experi- 
mental data, MC simulation with evaporation of light particles 
only and MC simulation including also the emission of IMF, 
respectively. 



(dashed histogram) which well reproduces the data of fig. 
I) and is characterized by an increased evaporation of IMF 
with respect to GEMINI. Again, to stress the importance 
of IMF emission, the dotted histogram shows the result 
of the simulation featuring the emission of light particles 
only. The reproduction of the experimental shape is still 
not perfect, but nevertheless the comparison shows the 
need for large perturbations of the kinematics, which in 
the present case have been obtained through the simula- 
tion of statistical emission of IMF. 



3.1.4 Correlation between AA and A 

In order to study the partition of the excitation energy 
between the two outgoing heavy fragments of the reac- 
tion, the experimental data were sampled in bins of re- 
constructed primary mass A of the PLF, for various win- 
dows of TKEL (corrected for Q gg , the Q- value between 
ground-states in the entrance and exit channels fl(i|| ). The 
ccntroids of the corresponding distributions of evaporated 
mass AA = A — A soc were then determined. The open 
squares and open circles in fig. [l(] present AA (without 
corrections) as a function of the primary mass of the PLF 
in the direct and reverse reaction, respectively. The data 
are shown for three bins of TKEL corresponding to partly 
damped events, where PLF can be safely distinguished 
from TLF due to the strongly anisotropic angular distri- 
butions jlj. The full symbols in fig. |l0| show the same 
experimental data after correction via Monte Carlo simu- 
lation. 

It has to be noted at this point that the applied correc- 
tions are indeed small (in most cases less than 1-2 amu). 
Moreover they are largely independent (within errors) of 
physical hypothesis (e.g. on energy partition), as it was 
checked by repeating the Monte Carlo simulations with 
different physical models. The corrections are largest (up 
to 4-5 amu) for the lightest masses at small TKEL, where 
the uncorrected data display, for the lowest masses, a kind 
of upward turn. This shape is the result of finite-resolution 
effects, mainly in the time-of-flight observable. It was al- 
ready pointed out by the authors of ref. pl[ that their 
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Fig. 10. Experimental average number of evaporated nucleons 
AA as a function of the primary mass A of the PLF in the 
direct (open circles) and reverse reaction (open squares) at 25 
AMeV, for three windows of TKEL. The full symbols show 
the experimental data after correction for the response of the 
setup, finite resolution effects and distortions of the analysis. 



eq. (8), which had been worked out for the analysis of the 
reconstructed primary mass distributions in the collision 
56 Fe + 165 Ho at 9 AMeV, just corrects each bin of a non- 
uniform distribution for the unequal contributions coming 
from the neighboring bins located to its right and to its 
left. Such a simple correction is justified only when the 
correlations among variables are negligible. In more com- 
plex cases, proper corrections can be performed only with 
the help of realistic Monte Carlo simulations, and this is 
the method applied in our analysis. 

Monte Carlo simulations showed that in our experi- 
ment, where the nuclei collide with large relative veloci- 
ties, at low values of TKEL the time-of-flight resolution of 
the PLF is the most critical parameter. Not only it mainly 
determined the resolution of the three variables used in 
fig. [lo] — namely A, A scc (used to build AA = A — A sec ) 
and TKEL — but, being a common ingredient of them 
all, it also introduced correlations among them. For ex- 
ample, A was positively correlated with A scc (leading to 
a partial compensation in AA) and with TKEL, however 
the details of these effects could be investigated only with 
Monte Carlo methods. 

The data of fig. |l0| present two distinct correlations be- 
tween the average number of emitted nucleons AA and the 
primary mass A for the two kinematic cases. When prop- 
erly corrected, these correlations appear to be approxi- 
mately linear and almost parallel. This finding is similar 
to the one already observed in the system 100 Mo + 120 Sn 
at 14 AMeV Q. Also in the present case, the most strik- 
ing feature resides in the different observed values of A A, 
depending whether they refer to PLF produced in the di- 
rect Nb+Sn or in the reverse Sn+Nb reaction. This is a 
direct evidence for a dependence of the average number of 
emitted nucleons on the net mass transfer. 
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3.2 Experimental results on light charged particles 

The purpose of the scintillator array "Le Mur" was to de- 
termine the average multiplicity of light charged particles 
(LCP) to be attributed to the PLF. As already stated, a 
first selection removed all particles stopped in the scintil- 
lator material. The remaining particles were then cleanly 
identified by charge Z, so that in the following we will 
refer just to hydrogen and helium ions, as their isotopic 
composition could not be determined with the employed 
experimental set-up. Due to the limited range of labora- 
tory angles (0i a b < 18.5°) covered with "Le Mur" in the 
present configuration, only a part of the angular distribu- 
tion of the light particles in the frame of the emitter could 
be measured. 

A first general selection of the data was performed by 
requiring that the light particles be emitted on the other 
side of the beam with respect to the PLF Silicon detectors, 
to avoid the complicated corrections due to their shadows. 
Further, in order to study the shape of the LCP velocity 
spectra, a forward angular range in the PLF frame was 
selected, such that the geometric acceptance of the scin- 
tillators did not appreciably bias the experimental data. 
This range was determined with the help of Monte Carlo 
calculations employing realistic velocity distributions, by 
requiring that the velocity distributions of the LCP be- 
fore and after the experimental filter had the same shape, 
within errors. It was found that the best angular range was 
#=18°-34° and 42°-54° for hydrogen and helium particles, 
respectively. 

In fig. O some experimental velocity spectra of LCP in 
the so determined angular range in the PLF frame (circles) 
are compared with the unfiltered results of evaporation 
calculations, obtained from the statistical code GEMINI 
(histograms) with appropriate excitation energies. The left 
and right columns of the figure are for hydrogen and he- 
lium ions, respectively. The first and second row refer to 
two windows of TKEL for the exit channel without net 
mass transfer in the collision 116 Sn + 93 Nb. The third 
row is for the exit channel leading to symmetry, either by 
a net loss of nucleons (left panel, mass ,4=104 in the col- 
lision 116 Sn + 93 Nb) or by a net gain of nucleons (right 
panel, mass .4=104 in the collision 93 Nb + 116 Sn). It can 
be seen that the agreement between experimental data 
and calculations is rather good thus showing that the es- 
timated excitation energies are correct and that, in the 
chosen angular range, the data are indeed rather unbiased 
and compatible with a thermal emission. 

Some trends, which will be explained in more detail in 
the following, are visible already in fig. O. The comparison 
of left and right columns shows that the velocity distribu- 
tions for hydrogen are generally broader than those for 
helium, having larger mean value, variance and skewness. 
The widths slightly increase with increasing TKEL (com- 
pare the first and the second row), whereas for a given 
TKEL they tend to increase or decrease for a net gain or 
loss of nucleons (compare the second and third row). Un- 
fortunately the extraction of an (apparent) temperature 
of the emitting source would be ambiguous. In fact, the 
AE vs. ToF identification technique does not give any in- 



formation about the mass of the detected particles. This 
fact prevents a reliable transformation of the velocity dis- 
tributions into energy distributions, especially for Z=l 
particles, where comparable contributions from protons, 
deuterons and tritons are expected. 

Now an analysis similar to that of sect. 3.1.4| can be 
performed also with the data of the light charged parti- 
cles Q. To do this, the average multiplicity (M\ cp ) of 
hydrogen and helium ions was determined as a function 
of the primary mass A of the detected PLF, for differ- 
ent windows on TKEL. In order to improve the statistics, 
a wider range from about 69=14° to 70° in the emitter 
frame was used. Monte Carlo simulations showed that set- 
up was still reasonably efficient in this range. In this way 
the relative multiplicity of light charged particles emitted 
by PLF of different primary masses A is relatively free 
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Fig. 11. Selected experimental spectra of the velocity of light 
charged particles in the PLF frame (circles) and results of evap- 
oration calculations with the statistical code GEMINI (his- 
tograms) . The left and right columns show results for hydrogen 
and helium ions, respectively. The first and second row are for 
two windows of TKEL for the exit channel without net mass 
transfer in the collision 116 Sn + 93 Nb, while the third row refers 
to the exit channel leading to symmetry, either by a net loss 
of nucleons (left panel, mass ,4=104 in the collision 116 Sn + 
93 Nb) or by a net gain of nucleons (right panel, mass A— 104 in 
the collision 93 Nb + 116 Sn). The widths a of the distributions 
were deduced from their second moments. 
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of bias, whereas the absolute values are somewhat more 
uncertain, because the extrapolation to the whole solid 
angle relies on assumptions about their in-plane and out- 
of-plane distributions. 

The left and right columns of fig. [l^ present the ob- 
tained results (efficiency corrected with Monte Carlo sim- 
ulations) for hydrogen and helium ions, respectively, in 
three windows of TKEL. The circles and squares refer to 
the multiplicities of light charged particles emitted from 
the PLF in the 93 Nb + 116 Sn and 116 Sn + 93 Nb reaction, 
respectively. The full symbols show the experimental data 
corrected under the assumption of a non-equilibrium ex- 
citation energy sharing (dependent on the net mass trans- 
fer), as deduced from the analysis of fig. [Io| (sec sect. L2 for 
details). The open symbols show the same experimental 
data after correction with an energy sharing mechanism 
independent of the net mass transfer. It is apparent that 
the results are very insensitive to the particular physical 
hypothesis on energy sharing used for the correction. 

The light charged particle multiplicity (Mi cp ) presents 
a general tendency to increase with increasing primary 
mass A of the emitting PLF. We want to focus attention 
on the fact that, in a given TKEL window, the multiplic- 
ities for the symmetric exit channel are quite different in 
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Fig. 12. Efficiency corrected experimental multiplicities for 
hydrogen (left column) and helium particles (right column) 
emitted by the PLF, plotted as a function of the primary 
mass of the PLF and for three windows of TKEL. The circles 
and squares refer to the multiplicities of light charged parti- 
cles emitted from the PLF in the 93 Nb + 116 Sn and 116 Sn + 
93 Nb reaction, respectively. The full and open symbols show 
the experimental data corrected assuming a non-equilibrium 
excitation energy sharing as deduced from the analysis of the 
data, and an energy sharing mechanism independent of the net 
mass transfer, respectively. The lines are drawn just to guide 
the eye. 



the two kinematic cases, although the mass of the emitting 
PLF is the same (A 104). They are larger in the direct 
reaction, where the PLF has gained mass, with respect to 
the reverse reaction, where it has lost an equal amount of 
nucleons. This finding is qualitatively in good agreement 
with the results already presented in fig. ^ for the total 
number of evaporated nucleons, AA, deduced on the base 
of a kinematic reconstruction. Thus also the emission of 
light charged particles shows a clear evidence for a corre- 
lation with the number of net transferred nucleons. This 
behavior is evident (beyond errors and uncertainties) at 
all TKEL values in case of helium ions, whereas for hy- 
drogens it becomes weaker with increasing TKEL. An in- 
terpretation of this fact will be given later in the paper 
(see sect. 4.3). 

It is finally worth noting that the measurements of 
light charged particles in the scintillator array "Le Mur" 
and of heavy reaction products in the gas detectors are 
completely independent. Therefore (M\ cp ) and A in fig. 
fl2| are free from the "instrumental" correlations affecting 



L4 and A, as discussed in sect. 3.1.4 



4 Discussion 

4.1 Energy sharing without net mass transfer 

A first information on the average partition of excitation 
energy between the two reaction partners can be deduced 
— in a substantially model-independent way — from the 
number of nucleons emitted in case of no net mass transfer, 
i. e. in the present work for A=93 (116) in the direct (re- 
verse) kinematics, respectively. The comparison between 
the two kinematic cases can be performed Q on the basis 
of a dimensionless parameter, Cf , representing the asym- 
metry in their global evaporation of nucleons FJ: 



C F = (AA h 116 - AA l g3 )/(AA h 116 



Ma) (1) 



where AA l g3 (AAi 16 ) is the total number of nucleons evap- 
orated from nuclei of primary mass A=93 (116) originat- 
ing from the entrance channel light (heavy) nucleus. As- 
suming a common value e for the average energy necessary 
to evaporate a single nucleon from the two reaction part- 
ners, Cf is also an estimate of the asymmetry in excitation 
energy partition (E^ 6 -E* l 3 )/(E^ 6 + E* l 3 ). The open cir- 
cles in fig. [b] show, as a function of TKEL, the value of 
Cf obtained from the present experimental data. In the 
same figure, the dashed line at Cf =0 and the dotted one 
at C F ~0.11 indicate the values expected for equal energy 
and equal temperature sharing, respectively. 

One can see that the so estimated average energy par- 
tition is fairly constant up to the highest accessible de- 
gree of inelasticity (TKEL « 700 MeV). If taken at its 
face value, the displacement of from zero (with the 
Sn fragments taking more than half of the total available 
excitation energy) would point to a situation intermedi- 
ate between equal energy and equal temperature sharing, 

The subscript F indicates a slightly different definition with 
respect to the parameter C of ref. |la| 
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Fig. 13. Experimental asymmetry Cf of the total mass evap- 
orated by the two reaction partners in case of no net mass 
transfer, as a function of TKEL (see text). The data refer to 
the asymmetric collision 93 Nb + 116 Sn at 25 AMeV. The open 
points show the data corrected with Monte Carlo calculations 
taking into account the effects of the analysis method. The 
full points represent the same data further corrected with the 
results obtained from an analysis of the symmetric systems. 
The dashed line at Cf — corresponds to the equipartition of 
excitation (equal energy sharing), while the dotted line shows 
the value expected for thermal equilibrium (equal temperature 
sharing) . 



even at low TKEL values. However, before drawing any 
conclusion, one should investigate the effects of a possible 
small dependence of e on the mass (and charge) of the de- 
caying nucleus. To do this, we took advantage of the data 
obtained, during the same experiment, for the two mass- 
symmetric systems 93 Nb + 93 Nb and 116 Sn + 116 Sn at the 
same bombarding energy of 25 AMeV. In fact, in the sym- 
metric exit channel of symmetric colliding systems, the 
PLF must take on average just half of the total excitation 
energy. Then one can readily estimate the average energy 
€93 spent by the 93 Nb nuclei to evaporate one single nu- 
cleoli, £93 = | TKEL/AA^ 3 m (and in a similar way en6 
for the 116 Sn nuclei). The superscript "sym" stresses that 
we refer here to symmetric colliding systems. The average 
cost to evaporate nucleons from Nb nuclei was £93 ~ 11-12 
MeV depending on excitation energy, while the value for 
Sn nuclei was found to be systematically lower by about 
0.7 MeV. All these values are compatible with statistical 
model calculations. 

It is perhaps worth noting that pre-equilibrium emis- 
sion, if present, might somewhat alter the absolute values 
of the so obtained e 93 and e ne, but to a much lesser extent 
their difference. Moreover, the effects of pre-equilibrium 
emission in the symmetric and asymmetric systems tend 
to cancel when applying such a correction to Cf- The so 
deduced small differences of e between Nb and Sn nuclei 
were used to correct the estimate of Cf in the asymmetric 
systems, by multiplying each term AA in cq. fy with the 
appropriate e value. The full circles in fig. [l^ show Cf vs. 
TKEL after correction. Within errors, the excitation en- 
ergy is now shared equally between the reaction products. 



In previous works in the literature p4| 17fl , referring 
to more asymmetric systems at lower bombarding ener- 
gies, the excitation energy was found to be almost equally 
shared between the colliding nuclei for small TKEL (cor- 



responding to a wide range of peripheral impact parame- 
ters) , while the thermal limit was slowly approached - but 
never reached - with increasing TKEL. Such a behavior 
has been qualitatively well understood in the framework 
of nuclear exchange models. Due to their randomness, the 
first exchanges of nucleons are predicted to produce about 
the same amount of excitation energy in the two cold nu- 
clei: the larger the mass asymmetry of the colliding nu- 
clei, the stronger the thermal non-equilibrium developing 
in this first phase of the interaction. However, this thermal 
disparity leads to a change of the fluxes of nucleons be- 
tween the partners of the reaction and drives the dinuclcar 
system towards thermal equilibrium. Therefore, allowing 
a long enough duration of the contact phase (as it may 
happen at small impact parameters and at low bombard- 
ing energies), the statistical decay of the two outgoing 
fragments should be consistent with thermal equilibrium 
(equal temperature condition). 

The data in this paper are the first on this subject 
at these higher bombarding energies. The striking feature 
is that Cf presents no evidence of an (even slow) trend 
towards the equal temperature limit (dotted line in fig. 
|l3|) with increasing TKEL. As such, the present data also 
differ from the results we obtained in a similar experiment 
at the lower bombarding energy of 14 AMeV (see fig. 2b 
of ref . H ) , where a trend towards the equal temperature 
partition was indeed found. 

An explanation can be sought in the strongly different 
interaction times rj nt at these two bombarding energies. 
As suggested by different model calculations at 25 AMeV 
the values of r; n t are shorter for a given TKEL and their 
increase with increasing TKEL is slower. Just as an ex- 
ample, the reaction times calculated in the frame of the 
Nucleon Exchange Model (lines) pcj and in a Landau- 
Vlasov approach (symbols) |3q| for the two bombarding 
energies are presented in fig. |14|. Referring, for example, 
to the NEM calculations, one sees that in the collision 
100 Mo + 120 Sn at 14 AMeV r int w 200 fm/c is obtained 
already at TKEL w 300-400 MeV, but similar interaction 
times are reached only at TKEL « 800-900 MeV in the 
collision 93 Nb + 116 Sn at 25 AMeV. The present indica- 
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Fig. 14. Reaction times calculated in the frame of the Nu- 
cleon Exchange Model (lines) and in a Landau- Vlasov ap- 
proach (symbols) for the systems 100 Mo + 120 Sn at 14.1 AMeV 
and 93 Nb + 116 Sn at 25 AMeV. 
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tion is in agreement with the result of ref . Q which claims 
that in the asymmetric system 40 Ar + nat Ag at 27 AMeV 
energy equilibration is not obtained for interaction times 
shorter than 180 fm/c. 

An experimental, although qualitative, observation in 
favor of a difference in time scales can be obtained also 
from the correlations d 2 a/dTKE ddp™ F (Wilczynski-Plot) 
measured at the two energies. At 14 AMeV, the high- 
est TKEL at which the parameter Cf could be analyzed 
corresponds to the attainment of an orbiting condition 
(implying long interaction times). In the experiments at 
25 AMeV, on the contrary, at the highest TKEL values 
used for determining Cf the PLF deflection angles are 
still peaked close to the grazing angle, indicating short 
interaction times even at TKEL as large as 800 MeV. 

Additional support to the hypothesis of a lack of sta- 
tistical equilibrium even at high TKEL can be obtained 
from the comparison of the data of the collisions 116 Sn + 
58 Ni, 116 Sn + 197 Au and 116 Sn + 116 Sn, all measured at 
25 AMeV. The strong mass asymmetry of the first two 
systems increases the sensitivity of the measurement of 
Z\Aii 6 on the possible evolution of the excitation energy 
sharing towards equilibrium. Let us consider, for all re- 
actions, the exit channel with Appp — WQ (no net mass 
transfer). Of course, in case of "equal energy" sharing and 
for a given TKEL, one expects the same excitation en- 
ergy of the PLF in the three systems. On the contrary, in 
case of an evolution towards statistical equilibrium ( "equal 
temperature" sharing), the excitation energy of the PLF 
is expected to depend on the system, as the mass ratios 
Aproj/A tot are very different (116/174«0.67 for Sn+Ni, 
116/313«0.37 for Sn+Au and 0.5 for Sn+Sn). Therefore, 
if AAuq depends just on the excitation energy of the PLF, 
it should increase differently as a function of TKEL. 

Figure 15 shows the experimental AAhq (uncorrected) 



for the two strongly asymmetric systems (squares and tri- 



¥40 - 

f 

3 

35 ; 

30 7 

25 ~ 7 

20 ' 7 

15 7 

10 7 

5 '- 



116 Sn+ 116 Sn 25AMeV 
□ 116 Sn+ 58 Ni 25AMeV 
116 Sn+ 197 Au 25AMeV 




100 200 300 



400 500 600 700 800 

TKEL (MeV) 



Fig. 15. Experimental evaporated mass AAne from excited 
PLF with A PLF = 116 for the three reactions 116 Sn + 58 Ni, 
93 Nb and 197 Au at 25 AMeV, as a function of TKEL. The full 
(dashed) line represents an empirical estimate of the AAne 
for an equal temperature sharing in the asymmetric Sn+Ni 
(Sn+Au) system, based on the data of the symmetric Sn+Sn 
system. 



angles) together with the results for the symmetric sys- 
tem 116 Sn+ 116 Sn (circles). The behavior expected for the 
asymmetric systems in case of "equal temperature" shar- 
ing (full and dashed lines for Sn+Ni and Sn+Au, respec- 
tively) can be obtained from the data points of the sym- 
metric Sn+Sn system by rescaling their abscissas with the 
ratios of total masses (116+197)/(116+116) and (116+58) 
/ (116+116), respectively. One sees that, up to the highest 
explored TKEL, the values of Z\An6 for the asymmetric 
system Sn+Au closely follow the points for the Sn+Sn 
system, thus suggesting an "equal energy" sharing behav- 
ior. Indeed, in the considered TKEL range, NEM calcula- 
tions give interaction times which are rather short (< 150 
fm/c) and very similar for the two systems. The behavior 
of the Sn+Ni system, on the contrary, is somewhat differ- 
ent and not fully understood, with data points in between 
the "equal energy" and the "equal temperature" expecta- 
tions. 



4.2 Energy sharing with net mass transfer 

Information on the dependence of the excitation energy 
sharing on the net mass transfer is carried by the slope of 
- and by the separation between - the two experimental 
correlations of fig. [Io| 

One easily sees that at 25 AMeV the difference of AA 
for a given A in the two kinematic cases amounts to about 
11-14 amu, almost independently of TKEL. If there are 
no violent dynamical effects at work, the emission of par- 
ticles is essentially of statistical nature and this difference 
can be ascribed mainly to a different excitation energy of 
the emitting PLF. This value of A A can be converted into 
an excitation energy using the previously found value of 
11-12 MeV for the average energy needed to remove one 
nucleon. Thus, for TKEL > 200 MeV (where the mass 
distributions of PLF in direct and reverse kinematics have 
some overlap), the observed difference suggests that a PLF 
of a given mass A produced in the direct Nb+Sn reaction 
has about 125 — 160 MeV of excitation energy more than a 
PLF of the same mass produced in the reverse Sn+Nb col- 
lision. In other words, there is an excitation-energy excess 
of about 6 MeV per net gained nucleon, averaged over the 
whole sequence of exchanges leading to the observed final 
TKEL. We recall that none of the usual ways of modeling 
the excitation energy sharing - neither the equal-energy, 
nor the equal-temperature scenarios, nor any combination 
of the two - foresees the observed existence in the corre- 
lation AA vs. A of two well-separated lines. 

For a given window of TKEL, the experimental corre- 
lations of fig. ^ can be treated, to a good approximation, 
as straight and parallel lines. This implies that the total 
number of nucleons emitted altogether, AA to t = AAppp + 
Z\AtlFi by any pair of reaction partners, Ap^p + ^4tlf = 
^4beam + ^target = ^tot , does not substantially depend on 
the particular exit channel. As in a previous work Q , the 
dependence of AA on the net mass transfer can be de- 
scribed by means of a dimcnsionless parameter, Rp , rep- 
resenting the asymmetry in evaporated mass for the exit 
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channel leading to symmetric division of the whole sys- 
tem 0: 



Ri 



(AA ; 
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sym 



sym m 



+ AA% m ) 



(2) 



where AA l sym (AA^ ym ) is the total mass evaporated from 
nuclei produced in the symmetric exit channel (that with 
primary masses A syul = Atot/2), originating from the light 
(heavy) colliding nucleus - and measured in our experi- 
ment as PLF in direct (reverse) kinematics. As already 
pointed out Rp is also an estimate of the excitation- 
energy asymmetry (E* ym - E*^) / (E* y]n + £*^ m ), in the 
limit that the small variations of e with excitation en- 
ergy can be neglected. Since the central mass A sym lies in 
the wings of the mass distributions where the statistics is 
lower, Rp was not obtained directly from the data, but 
deduced from the result of a simultaneous linear fit to the 
parallel correlations of fig. |lG. 

Using the expressions cq. |l| and || for Cf and Rf, and 
with a common average value of e, the excitation energy 
for products of primary mass A 1 (A h ) deriving from the 
original light (heavy) colliding nucleus can be cast in the 
form §]: 
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where A l (Aq) is the lighter (heavier) mass between A oe 



and A taTgct in the entrance channel 
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experimental slope includes a contribution (that with the 
term Cf) simply describing the dependence of excitation 
energy on mass, while only the term with Rp truly repre- 
sents a dependence on net mass transfer and, as such, it is 
responsible for the existence of two distinct correlations. 

The circles in fig. |l6| a show the obtained values of 
i?F/Aiif as a function of TKEL, together with the results 
for the systems 100 Mo + 120 Sn at 14 AMeV (triangles) 
and 74 Ge + 165 Ho at 8.5 AMeV (crosses). In this latter 
case, where only the correlation for the direct reaction had 
been measured, the points were reconstructed on the ba- 
sis of the measured slopes and of the average partition as 
estimated in ref. []l5[). In all three cases the percentage of 
excitation energy gained by the net transfer of one nucleon 
is strongly decreasing with TKEL. 

If we consider the product (Rp /A^f) • TKEL, which 
is the coefficient of the net-mass-transfer dependent term 
in eq. ^, we obtain the result shown in fig. |l6|b. At 14 and 
25 AMeV this product is roughly constantfor it slightly 
increases with TKEL), typical values at TKEL<550 MeV 
being 2.5-3 and 5.5-6.5 MeV per net transferred nucleon, 
respectively. Although only two bombarding energies may 
be insufficient to draw any stringent conclusion and the 
observed factor of about 2 between the two experiments 
might be rather fortuitous, it is worth noting that they 

2 The subscript F stresses the different definition with re- 
spect to the parameter R of ref. [Jl5| 



scale like v^ el at contact. This observation agrees with the 
expectations of exchange models where the recoil momen- 
tum of the transferred nucleon(s) is the main contribution 
to the energy dissipated in a single exchange process and 
the observed correlation with the net mass transfer may 
arise from an intrinsic asymmetry in the excitation en- 
ergy generated in the donor and acceptor nucleus at each 
elementary step. Also BNV calculations |37| at 14 AMeV 
ascribed such a correlation between total evaporated mass 
and net mass transfer to an intrinsic asymmetry in the nu- 
cleon exchange process. 

It has to be noted that in the frame of an exchange 
picture, one would also expect a flattening |3£| of the 
slope of the correlations AA vs A with increasing TKEL 
(as observed indeed in the collision 74 Ge + 165 Ho) and 
hence a decrease of the coefficient (Rp /A^n) ■ TKEL. In 
fact, for a given net mass transfer, with increasing TKEL 
there should be a growing contribution from exchanges 
taking place at later stages of the reaction, when the rel- 
ative motion is somewhat slowed down (and hence the 
dissipated energy per exchanged nucleon is lower). How- 
ever, rather surprisingly, the experimental data show that 
(i?F/^4dif) - TKEL does not decrease, on the contrary it is 
constant or even a weakly increasing function of TKEL. 

In the frame of an exchange picture, one would like to 
have access to the asymmetry in excitation energy at each 
elementary step: 
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Fig. 16. a): Percentage of excitation energy gained by the 
net transfer of one nucleon, Rp/Ada, for the systems 93 Nb + 
116 Sn at 25 AMeV (circles), 100 Mo + 120 Sn at 14.1 AMeV (tri- 
angles) and 74 Ge + 165 Ho at 8.5 AMeV (crosses), as a func- 
tion of TKEL. b): Coefficient (R F /A dit ) ■ TKEL of the net- 
mass-transfer dependent term in the excitation-energy versus 
primary-mass parametrization (see text). Same symbols as in 
part a). 
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where e a and are the excitation energy generated in 
the acceptor and donor nucleus, respectively. Based on 
some assumptions to be discussed later, Toke at al. jl5| 
proposed an original method to estimate ?y from the ex- 
perimental data. Their prescription (see eq. A6 of p5[) 
can be expressed, employing our parameter notation, in 
the following way: 



Adif 



£- al TKEL) 



?7exp — 



7(2) 



2 , TKEL) 



C F {2)-C F {\) 



(5) 



where a \ is the variance of the mass distribution and the 
indices (1) and (2) refer to the experimental data taken 
from two successive TKEL bins, with TKEL(2)>TKEL(1). 
The quantity ?7 CX p is thus, by construction, a differential 
quantity pertaining to a particular stage of the collision. 

In the case of 74 Ge + 165 Ho at 8.5 AMeV, Toke et al. 
found that their prescription ry exp gave a reasonable value 
of « 0.3 for the first 100 MeV of TKEL, increasing to 
about 1 at 130 MeV. This corresponds to an asymmet- 
ric sharing of excitation energy in the ratio 1:2 between 
the donor and the acceptor nucleus. The results of the 
same prescription applied to the systems 93 Nb + 116 Sn at 
25 AMeV (circles) and 100 Mo + 120 Sn at 14 AMeV (tri- 
angles) are displayed in fig. [l7] together with the data for 
the original 74 Ge + 165 Ho system at 8.5 AMeV (crosses). 
In our systems the obtained 77 exp starts at values around 1 
for the lowest TKEL and systematically increases with in- 
creasing TKEL, up to values of about 10-20 for the largest 
TKEL. This fact, of course, prevents the interpretation of 
?7exp as an estimate of the excitation energy asymmetry in 
a single exchange process. Such a quantity is indeed ex- 
pected to be bound between +1 and -1, unless one accepts 
the (unconventional) hypothesis that one of the two nu- 
clei, acceptor or donor, experiences on average a negative 
variation of excitation energy (i.e., the nucleus becomes 
systematically colder in an exchange process). 
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Fig. 17. Value of the experimental quantity Jfcxp, which, in 
a stochastic nucleon exchange picture, is expected to estimate 
the excitation energy asymmetry in an elementary exchange 
process, as a function of TKEL, for the systems 93 Nb + 116 Sn 
at 25 AMeV (circles), 100 Mo + 120 Sn at 14.1 AMeV (triangles) 
and 74 Ge + 165 Ho at 8.5 AMeV (crosses). 



One is then lead to reexamine critically the assump- 
tions on which the derivation of eq. ^ relies, due to the 
fact that at high TKEL r\ appears to be too large by 
about a factor of 10. Basing on the one-body exchange 
picture (which portrays the collision as an evolutionary 
process consisting of a long sequence of elementary steps, 
each producing a small change of the system) the main 
assumptions are: 

i) " collisions leading to different values of i?i oss 

evolve along a common path in the space of the relevant 
collision parameters..." |lq| ; 

ii) in non-central collisions the dissipation mechanism 
proceeds mainly via the exchange mechanism, with other 
inelastic or collective excitation modes playing a minor 
role, if any; 

iii) the dependence of the excitation energy sharing on 
net mass transfer is quite well approximated by a linear 
relationship, like for example that of eq. |3|; 

iv) the exchange mechanism proceeds through stochas- 
tic transfer of single uncorrelated nucleons and the num- 
ber of exchanges can be deduced from the variances of 
the mass distribution, V cxc h = o\. 

Point i) implies that the evolution of the relevant quan- 
tities (namely, for the present discussion, the number of 
exchanges, the mass variances and the excitation energy 
sharing parameters) proceeds along a common path for 
collisions leading to consecutive bins of TKEL. In other 
words, the evolution of collisions leading to TKEL and 
TKEL+ZiTKEL differs only in the final part of the pro- 
cess, the one which is responsible for the additional dis- 
sipation ZiTKEL. This is of course an approximation, as 
"...different values of -Eioss are expected to involve differ- 
ent partial waves in the entrance channel and, therefore, 
different system trajectories" jlj|. In particular, in the 
derivation of eq. ^| it is assumed that the number of ex- 
changes in the evolution from TKEL(l) to TKEL(2) is 
simply given hy a\(2) ~ a\{\) . 

In order to check how far this approximation is justi- 
fied, the code for the Nucleon Exchange Model by Ran- 
drup [ p0| was run with the following minor modification. 
In the original code, the equations of motion of several 
relevant collective variables are numerically integrated in 
short time steps along the trajectory of the system and 
only at the end of the interaction the final values of mass 
(and charge) variances are estimated. The applied modifi- 
cation performs an estimation of the current mass variance 
at each integration step, so that it is possible to follow the 
evolution of u\ along the (average) trajectory leading to 
any given final dissipation TKEL. In particular, it is then 
possible to judge in how far the mass variance cr\(l) of a 
trajectory with final dissipation TKEL(l) is a good esti- 
mate of <j\(2 |TKEL(1)), that is the mass variance of the 
trajectory with final dissipation TKEL(2) at the moment 
when the dissipation had reached the intermediate value 
of TKEL(l). This is shown in fig. [l8], where each curve 
labeled with the value of the final TKEL in exit channel 
displays the evolution of a\ with the dissipated energy 
along the trajectory. If the hypothesis of i) were rigor- 
ously correct, all these curves should fall one on top of the 
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Fig. 18. Evolution of a\ with the dissipated energy along se- 
lected trajectories leading to the indicated final TKEL values, 
calculated with the Nucleon Exchange Model. 



other. One sees that when estimating a\ in a previous in- 
stant on a definite trajectory by means of the final mass 
variance of a lower curve, one introduces a systematic er- 
ror of the order of 10% on o\{\). This corresponds to an 
overestimation of 77 eX p of about 20-30%, at most, and it is 
therefore insufficient to explain by itself the surprisingly 
high value of r^xp- 

The importance of inelastic excitations in the two col- 
liding nuclei (point ii) — either of single particle type or 
of more collective nature, but not strictly related with a 
mechanism of nucleon exchanges — has been assessed by 
some authors in the literature especially for the first phase 
of the reaction. Assuming that the energy dissipation pro- 
ceeds partly via such an exchange mechanism and partly 
via inelastic excitations (not dependent on the net mass 
transfer), one can still obtain an experimental value of 
?7exp ■ However in this case only part of the average increase 
of excitation energy of the nucleus, (E*(2)) — {E*(l)} cx 
Cp(2) — Cf(1), is due to the exchange process (this corre- 
sponds to addinga term -E* nelastic to the right hand side of 
eq. (A5) in ref. p5j)- Thus the values of r^xp presented in 
fig. O, where all dissipation was attributed to an exchange 
mechanism, would be a lower limit of the excitation energy 
asymmetry of a single exchange rj. 

The hypothesis iii) of an approximately linear depen- 
dence of the excitation energy sharing on the net-mass- 
transfer seems quite well supported by the present data 
(see fig. [lO]) as well as by the data of the collisions 100 Mo 
+ 120 Sn at 14 AMeV § and 74 Ge + 165 Ho at 8.5 AMeV @. 

It is commonly assumed (hypothesis iv) that the vari- 
ance a\ of the experimental mass distributions is a good 
estimator of the number of elementary exchanges, an as- 
sumption which is justified in the frame of a mechanism 
of stochastic exchanges of single nucleons. It is well known 
that at low energy losses the results of this procedure are 
in rather good agreement with the number of exchanges 
predicted by theoretical calculations. However, with in- 
creasing TKEL, there is a growing discrepancy, as the 
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Fig. 19. Comparison of the experimental mass variances (sym- 
bols) with those calculated with the Nucleon Exchange Model 
(lines) for the systems 74 Ge + 165 Ho at 8.5 AMeV §g] (aster- 
isks), 100 Mo + 120 Sn at 14.1 AMeV (triangles) g] and 93 Nb + 
116 Sn at 25 AMeV (circles). 



mass variances tend to increase more rapidly than the 
theoretically calculated number of exchanges. For exam- 
ple, in the system 74 Ge + 165 Ho at 8.5 AMeV §§, al- 
ready at TKEL= 100 - 150 MeV the experimental mass 
variances become about 2-3 times larger than predicted 
by NEM. In our systems, where much larger amounts of 
energy are dissipated, the discrepancy between experimen- 
tal and calculated mass variances becomes dramatic, as it 
amounts to about a factor of 5 at the highest TKEL values 
where r) exp attains values of about 10. The experimental 
variances in comparison with NEM calculations for these 
three systems are shown in fig. 

Using the mass variances calculated in the model would 
give values of rj cxp in a meaningful range, but the discrep- 
ancy with the experimental mass variances is huge and 
remains to be understood. One might try to explain the 
discrepancy by assuming that the flow directions of succes- 
sively exchanged nucleons are to some extent correlated, 
or even that clusters of nucleons (instead of single nucle- 
ons) are transferred in a single exchange process. In this 
latter extreme case, one should divide the mass variance 
by the average mass \i of the exchanged cluster, o\j \i, 
in order to have an estimator of the number of elemen- 
tary exchange steps. However, it seems difficult to accept 
values as large as /i= 4-5 amu, which are needed at the 
highest TKEL. 

In conclusion, remaining in the frame of an exchange 
picture it seems difficult to find a way to bring the high val- 
ues of the experimental quantity 77 exp down into a reason- 
able range, where it can be interpreted as a good estimator 
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of the excitation energy asymmetry in an elementary ex- 
change process, rj. Recalling that, conceptually, eq. |^ is 
worked out in the frame of a one-body exchange picture 
(and it looses meaning outside of that frame), this failure 
casts doubts on the validity (or the relevance) of such an 
exchange picture for describing the dissipation mechanism 
at these higher bombarding energies. 

In any case, one should keep in mind that also other 
mechanisms, like fluctuations in the neck rupture, may 
contribute in a nonnegligible way to the excitation energy 
sharing and to the width broadening of the experimental 
mass distributions. 



4.3 Angular momentum sharing and light charged 
particles 

The ex peri mental data on light charged particles presented 
in sect. 3.2 show evidence for a correlation with the num- 

similar to that of the total 



A 
I 



ber of net transferred nucleons, 
evaporation from the heavy fragments. 

The multiplicity of light charged particles is an increas- 
ing function of the excitation energy of the emitting nu- 
cleus. However, due to the limited solid angle covered in 
the present experiment by the detector array "Le Mur" , 
the absolute values of the light charged particle multiplic- 
ities are affected by uncertainties larger than those associ- 
ated with their ratios. Thus our analysis has been focused 
on the ratio between the average multiplicities of Hydro- 
gen and Helium particles, (Mn) / (Mn e ) , emitted from the 
detected PLF 

The experimental data for particles emitted by the 
PLF in the direct (full circles) and reverse (full squares) 
collision of the system 93 Nb + 116 Sn at 25 AMeV are pre- 
sented in fig. pO| a for exit channels without net mass trans- 
fer. In agreement with the results on th e average excitation 
energy partition presented in sect. LI (see also fig. [l^), an 
equal division of the total excitation energy (estimated by 
TKEL) has been assumed, although the arguments that 
follow are rather insensitive to this hypothesis. The two 
sets of experimental data are very similar, possibly be- 
cause of a weak dependence of the multiplicities of light 
charged particles on the mass of the emitting nucleus in 
this mass region. 

On the same figure, the results of evaporation calcu- 
lations with the statistical code GEMINI Q are shown by 
the open symbols for a 116 Sn nucleus in case of zero spin 
and large spin. In this latter case, the calculations show 
also the additional effect of a prolate deformation (with 
a representative axis ratio of 1.6), which for large spins 
may be more appropriate than a spherical shape. The ra- 
tio between Hydrogen and Helium particles appears to be 
sensitive to the angular momentum of the evaporating nu- 
cleus, with large angular momenta (and consequent defor- 
mation) favoring the emission of the more massive Helium 

3 The latest version (August 2000) of the code has been used, 
with Hauser-Feshbach formalism for light particles up to Li 
(Z_imf_min=4), including the ground-unstable particles 5 He 
and Li (exotic_index=2), and transmission coefficients from 
the incoming wave boundary condition model (tl_iwbc=.true.) 
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Fig. 20. a): Ratio of the experimental average multiplicities 
of Hydrogens and Helium particles, (Aiu) / (Mhc) , emitted by 
the PLF in exit channels without net mass transfer. The full 
circles (squares) refer to PLF from the direct (reverse) collision 
93 Nb + 116 Sn at 25 AMeV and are plotted as a func tion of the 



excitation energy estimated from the data (see sect. 4.3 in the 
text). For comparison also the results of Gemini calculations 
for two spin values of an emitting spherical 116 Sn-source are 
shown. In case of high spin, the calculations were done also for 
a prolate deformed source, with an axis ratio of 1.6. b): Same 
presentation as in part a) except that the experimental data 
refer to events leading to a symmetric mass division in the exit 
channel and the calculations to a 105 Pd-source. 



particles with respect to the lighter Hydrogens. Thus, the 
experimental rapid drop of (Mn)/{Mn e ) with increasing 
E* can be ascribed to the rise of the average angular mo- 
mentum of the emitting nucleus. Indeed, the fraction of 
orbital angular momentum transferred into spin of the col- 
liding nuclei is small in peripheral collisions and increases 
when going to more central collisions (that is to larger 
TKEL values) 

Figure po| b presents the ratio of light charged particles 
{Mn)/{Mn e ) emitted by the PLF for events leading to 
symmetric mass division in the exit channel. Again, the 
two sets of experimental data refer to PLF measured in 
the direct and reverse kinematics for the collision 93 Nb + 
116 Sn at 25 AMeV (full circles and squares, respectively). 
In this case, the excitation energy E* of the PLF has been 
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estimated from the measured TKEL assuming an excita- 
tion energy division in agreement with th e find ings on 
the evaporated number of nucleons (see sect. 3.1.4 and fig. 
|Io| ). The horizontal bars show the range of excitation en- 
ergies spanned by energy divisions between the presently 
adopted one and that for equal energy sharing. For com- 
parison, GEMINI calculations for a nucleus of 105 Pd are 
also drawn. 

As already explained in rcf. pdfr , for PLF measured 
in the direct reaction — which therefore experienced a net 
mass gain of nucleons — the experimental results (circles) 
indicate rather low ratios (Mn) / {M-^e) , thus pointing to 
high spin values. The opposite holds for PLF measured 
in the reverse reaction — hence produced by a net loss of 
nucleons — where larger ratios (Mn)/{Mn e ) indicate lower 
spins of the emitters. The same conclusion would hold true 
also in case of different excitation energy partitions, indi- 
cated by the horizontal bars of fig. po|b. Thus, with respect 
to the expectations for full equilibrium, the net gain (loss) 
of nucleons seems to be correlated with an excess (reduc- 
tion) of both excitation energy and angular momentum 
sharing. 



5 Conclusions 

The collision 93 Nb + 116 Sn at 25 AMeV has been studied 
in direct and reverse kinematics. The analysis of primary 
and secondary masses of the PLF demonstrates the exis- 
tence of two distinct correlations between total evaporated 
mass AA and primary mass A in the two kinematic cases. 
This shows that the total evaporated mass depends on the 
net mass transfer between the two colliding nuclei. Also 
the data concerning the multiplicity of light charged par- 
ticles present a similar difference between the two kine- 
matic cases. Both experimental observations, which are 
independent of each other, can be interpreted in terms of 
a dependence of the excitation energy sharing on the net 
mass transfer. Moreover, the measured ratios of Hydrogen 
and Helium multiplicities strongly point to a correlation 
also of the angular momentum sharing with the net mass 
transfer. 

An explanation of the experimental findings in the 
frame of a stochastic exchange picture is challenged by 
the persisting strength of the correlation between AA and 
A even at high TKEL. The failure to obtain from the data 
a meaningful estimate of the excitation energy asymme- 
try r\ in an elementary exchange process casts doubts on 
the validity (or the relevance) of such an exchange picture 
for describing the dissipation mechanism at these higher 
bombarding energies. Indeed, with increasing TKEL and 
bombarding energy, other mechanisms may become im- 
portant, which cannot be described simply with the el- 
ementary process of exchanging matter across a window 
between the two nuclei. For example, a relevant role could 
be played by dynamic effects, such as formation and rup- 
ture of a neck during the collision, a mechanism which in 
principle might explain both the correlation of excitation 
energy with net mass transfer and the very large exper- 
imental widths of the mass distributions. A quantitative 



estimate of such collective effects is beyond the scope of 
the present work and requires extensive and detailed the- 
oretical calculations. 
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discussions. 



Appendix: Monte Carlo simulations 

Details of the general structure of Monte Carlo simula- 
tions were already given in sect. 2.3 of ref. [Q, while the 
Appendix A of ref. (To) was more concerned with the 3- 
(and 4-)body events and the typical resolutions obtained 
in their reconstruction. 

Briefly, a binary dissipative collision is simulated on 
the basis of realistic distributions for the relevant physical 
quantities. The Total Kinetic Energy (TKE) in the exit 
channel is considered as the leading variable to describe 
the reaction mechanism. In fact, according to theoreti- 
cal models and to the usually adopted picture (see e.g. 
fl4l|), one expects — on average — a monotonic correla- 
tion of TKE with both interaction time Ti n t and angular 
momentum t m in the entrance channel (TKE decreases 
with increasing Tj n t and with decreasing t ln ). Therefore, 
first of all a realistic distribution of TKE is generated and 
randomly sampled. 

Mean values and variances of the masses and scattering 
angles of the reaction products are then empirically pa- 
rameterized as a function of TKE and these parametriza- 
tions are iteratively tuned until realistic correlations - 
similar to the experimental ones — are produced, both 
for TKE-mass {diffusion-plot) and TKE-# cm {Wilczynski- 
plot). 

The value of TKE is used to estimate the total excita- 
tion energy of the system (after correcting for an average 
Q-value between entrance and exit channel). The total 
excitation energy can then be divided between the two 
reaction partners either equally {equal energy sharing) or 
in proportion to their masses {equal temperature sharing) 
or in agreement with the excitation energy division ex- 
perimentally deduced from the total number of nucleons 
evaporated from the fragments. 

The assumed monotonic decrease of TKE with l\ a (or 
equivalently with impact parameter 6), allows to extract 
also an average correspondence between TKE and £i n , via 
the relation f^°™ {da rcac /dE) AE = vrX 2 (^ r _ £fj w here 
£ gT is the grazing angular momentum corresponding to 
E cm . From theoretical models E^,^l|, an estimate of the 
dissipated angular momentum (that is of the amount of 
angular momentum transferred from the orbital motion 
into the internal degrees of freedom of the interacting sys- 
tem) can be obtained. The dissipated angular momentum 
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can then be divided between the two reaction partners 
cither proportionally to their moments of inertia (as ex- 
pected from a more equilibrated process) or with an excess 
in the nucleus gaining mass (as suggested by the present 
experimental results). 

At this point, all the primary quantities of the simu- 
lated binary dissipative collision are defined and the next 
step is to simulate the decay of each fragment. 

It was noted |lj that, for the purpose of studying just 
the heavy primary fragments, the particle emission in the 
Monte Carlo simulations needs to be nothing more than 
a way of adding statistical perturbations to the primary 
velocity vectors. As the results were not very sensitive to 
the adopted multiplicity distributions or relative yields of 
light particles, a simple rough paramctrization was then 
good enough and no big effort was devoted to tuning the 
evaporation step in a realistic way. 

However, in the present work, where not only the heavy 
reaction products but also the light charged particles are 
studied, it becomes important to perform a better simu- 
lation also of the evaporation step. 

For each decaying primary fragment, the dependence 
of the evaporation step on the initial values of excitation 
energy E* , spin J and mass A of the emitter is modelled 
according to statistical model calculations performed with 
the code GEMINI p7|. However, a direct event-by-event 
coupling of the Monte Carlo simulation code with GEM- 
INI results unpracticable because too time-consuming. In 
fact, filtered simulated data with statistics comparable 
or even larger than the experimental data are desirable 
(to keep additional statistical fluctuations in the results 
low), but then the required computing time becomes un- 
affordablc. Instead, a series of GEMINI calculations has 
been performed for a grid of values of E* , J and A of 
the emitter. The multiplicities of the various evaporated 
light particles and intermediate mass fragments have been 
acco rding ly parametrized. Actually, as already stated in 
sect. 3.1.31, it is necessary to tune the IMF multiplicities 



predicted by GEMINI in order to better reproduce the 
experimental spectra of AA and Zikcm shown in fig. |8|and 
^, respectively. At each step of the decay chain, the species 
to be emitted is randomly chosen according to the rela- 
tive importance of such multiplicities. Actually, in order to 
mimic the correlations among successive steps of the de- 
cay chain (which are naturally displayed by the full Gem- 
ini calculations), the random choice is further weighted so 
as to favor decay steps leading to daughter nuclei nearer 
to the so-called EAL line Q. 

At the moment of the emission (that is at the exit 
channel barrier), each particle or intermediate mass frag- 
ment has a kinetic energy of thermal origin. In order to 
speed up the simulation, this energy is sampled from a 
distribution with a (surface) Maxwellian shape, the ap- 
propriate temperature being that in the daughter nucleus 
at the barrier. Of course, to obtain the asymptotic kinetic 
energy of the particle it is necessary to add the additional 
contribution due to the Coulomb repulsion. 

For each fragment, the evaporation process is followed 
along all the decay chain until the excitation energy of 



the decaying nucleus is almost completely exhausted and 
the 7-ray decay becomes dominant (near the threshold for 
particle emission). 
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